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We present a quantitative analysis of an escape rate for switching from the superconducting 
state to a resistive one in series arrays of globally coupled Josephson junctions. A global coupling 
is provided by an external shunting impedance. Such an impedance can strongly suppress both 
the crossover temperature from the thermal fluctuation to quantum regimes, and the macroscopic 
quantum tunneling (MQT) in short Josephson junction series arrays. However, in large series arrays 
we obtain an enhancement of the crossover temperature, and a giant increase of the MQT escape 
rate. The effect is explained by excitation of a spatial-temporal charge instanton distributed over a 
whole structure. The model gives a possible explanation of recently published experimental results 
on an enhancement of the MQT in single crystals of high- Tc superconductors. 
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Great attention has been devoted to an experimental 
\ and theoretical study of do biased series arrays of Joseph- 
son junctions [E 0, IE 01 ■ Such a system displays diverse 
fascinated nonlinear classical and macroscopic quantum- 
mechanical phenomena. E.g. a resistive state of Joseph- 
son junction series arrays can show synchronized behav- 
ior Q, 0, and this effect has been used in Josephson 
\ voltage standard devices As we turn to a region of 
small dc bias currents and low temperatures, the macro- 
scopic quantum-mechanical phenomena start to play a 
role. Thus, a quantum phase superconductor-insulator 
transition has been observed in artificially prepared se- 
ries arrays of a small size Al / AI2O3 / Al junctions All 
these effects strongly depend on the interaction between 
Josephson junctions. 

This field of research, i.e. the macroscopic quantum 

[ phenomena in spatially extended superconducting sys- 
tems, has been boosted even further by recent discov- 

\ ery of macroscopic quantum tunneling (MQT) in single 
crystals of layered high-Tc superconductors H, 0- At 

' low temperatures the MQT determines the escape rate 
of the switching from the superconducting state to a re- 
sistive one. Although the MQT of a single Josephson 
phase has been found long time ago in low- Tc lumped 
Nb Josephson junctions 0,01, MQT in layered high-Tc 
superconductors has shown many novel features. A most 

I unexpected result is that the MQT escape rate Tmqt 
in high-Tc superconductors is four order of magnitude 
larger than the MQT escape rate for a lumped Joseph- 
son junction having the same parameters 0|. Moreover, 
the crossover temperature T* from the thermal fluctua- 
tion regime to the MQT regime increases in respect to 
a lumped Josephson junction. In these experiments it 
was also found that the escape rate Tt in the thermal 
fluctuation regime did not differ from the escape rate of 
a single Josephson junction. 

Layered high-Tc superconductors can be modelled as a 
stack (a series array) of intrinsic Josephson junctions 



A modern fabrication technique allows to prepare single 
crystals of layered high-Tc superconductors with an ex- 
tremely homogeneous distribution of critical currents of 
intrinsic Josephson junctions, and a low level of dissipa- 
tion . Since in the model of independent Josephson 
junctions the escape rate F is just proportional to A'^, 
an enhancement of the MQT observed in layered high- 
Tc superconductors stems from an interaction between 
intrinsic Josephson junctions. All experimental obser- 
vations receive a natural explanation in a simple model 
IQj of Josephson junctions series array with an intrin- 
sic charge interaction between nearest-neighbor Joseph- 
son junctions 0, 0|. Moreover, such a model allows 
quantitative comparison with experimental results, and 
a good agreement has been found as the Debye screening 
leng th is of the order of superconducting layer thickness 

EJ 

However, the authors of Ref. :6| proposed an other 
model in order to explain a giant increase of the MQT 
escape rate. In this model intrinsic Josephson junctions 
are globally coupled due to the presence of an external 
shunting impedance Z (see schematic in Fig. 1). An 
influence of electromagnetic environment, and in partic- 
ular, a shunting impedance Z on the MQT in a lumped 
Josephson junction has been studied long time ago in 
Refs. 0, Q, 01 . It was shown that the presence of a 
small shunting impedance can lead to a strong suppres- 
sion of the MQT in a lumped Josephson junction. There- 
fore, natural questions arise: what is a role of shunting 
impedance Z in the macroscopic quantum dynamics of 
large Josephson junction series arrays {N >> 1) and 
what is a most appropriate model in order to explain 
a giant increase of the MQT escape rate in Josephson 
junction series arrays? 

In order to answer these questions we carry out a quan- 
titative analysis of the escape rate F(/) in globally cou- 
pled Josephson junction series arrays. A global coupling 
is provided by an external shunting impedance Z (see Fig. 
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FIG. 1: Schematic of a dc biased layered high- Tc supercon- 
ductor and a series array of Josephson junctions. A strongly 
localized instanton (dashed line) and a charge instanton with 
long tails (solid lines) are shown. 



1). We obtain, and it is a main result of the paper, that 
if the electromagnetic environment strongly suppresses 
the MQT in a lumped Josephson junction, the role of 
an external shunting impedance is diminished in large 
Josephson junction arrays, and the standard quantum- 
mechanical behavior is recovered. Therefore, the glob- 
ally coupled Josephson junction series arrays can show a 
giant increase of the MQT escape rate with a strong de- 
pendence on a number of junctions. Moreover, the MQT 
escape rate is tunable in a wide region by a simple change 
of Z. Such tuning of the MQT escape rate can be very 
promising for a modern field of quantum information pro- 
cessing [l6( . 

A Josephson junction series array is characterized by 
the set of Josephson phases (^„(t), where number n 
changes from 1 to A^. Moreover, electrodynamics of a 
shunting impedance Z is described by a flowing charge 
Q{t). We will consider a particular case as all resistive 
effects are small, i.e. the shunting impedance has only a 
reactive part, and it contains inductor L and capacitor 
C in series. In order to obtain the escape rate F we use 
"instanton technique" 0,0,0], and therefore, (pn{T) 
and Q{t) arc periodic functions of the imaginary time r 
varying from to Ti/{kBT) (T is the temperature). The 
escape rate is determined by the action S as 



F ~ y" Dq{T)D^n{T)exp 



S{q,ip„} 



h/kBT 



L{T)dT 



(1) 



and the Lagrangian of a series array with the shunting 
impedance is written as 



V L/„ -I- i — V q{T)ip„ (r) 



Uni^p) = C0S((3„(t) + jipnir) , 



(2) 



where j = I/Ic is the normalized external dc current, 
and Ic is the nominal value of the critical current of a 
single junction. Here, Up is the plasma frequency of a 
single Josephson junction in the absence of dc bias. The 
Lagrangian is expressed in units of Ej^ where Ej is the 
Josephson energy of a single junction. The q = Q/^EjC 
is the normalized charge flowing through the impedance 
Z . The shunting impedance is characterized by the res- 
onance frequency = l/^J LC , where L and C are the 
impedance inductance and capacitance, accordingly. The 
coupling between the Josephson junction series array and 
the shunting impedance branch is described by parame- 
ter a — yjLj/ and Lj = 2e/{hlc) is the Josephson 
inductance. 

Integrating over we obtain the effective 

action Sefj that depends on the variables c^n('^) only 

h 
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where the kernel Gt{j) is determined as 



(3) 



Gt[t) 



knT 



E 
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m{2TTkBT)/h, m = ±l,±2.. 



(4) 



Thus, the last term in Eq. © presents an effective global 
charge interaction, that is due to current fluctuations 
flowing through an external shunting impedance. 

In the escape experiments Ej » hujp, and the switch- 
ing to a resistive state occurs as the dc current / is close 
to Ic, and therefore (j — 1) <C 1. In this case the potential 
Uni^p) is written as 



Un{ip) = (1 -j)(PniT) 
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(5) 



The escape rate is determined by the particular solution 
</j„(r) providing the extremum of effective action ((JJ. At 
high temperatures such a solution is determined by ex- 
tremum points of the potential J7„, and it is written as 



= 2^2(1 -j)5„, - v/2(l-j) . 

Here, I is a junction number where the fluctuation oc- 
curs. Since this solution does not depend on the time r, 
we can immediately conclude that the last term in (O 
does not give contribution to the escape rate exponent 
Ft — exp(— S'{(/3j}/?i. Note here that an absence of 
the dependence of the escape rate exponent in the ther- 
mal fluctuation regime on a number of junctions A^ is 



3 



a generic property of Josephson junctions series arrays 
with a charge interaction |l2|. 

However, the crossover temperature from the thermal 
fluctuation regime to the MQT can be strongly enhanced 
by such a global coupling. Indeed, using the method elab- 
orated in |l3 . lrH . ll8| we obtain that at high temperatures 



the optimal fluctuation 
has a form: 



(r) around an extremum point 



(6) 



where the eigenfunctions 0„ are the solution of the non- 
local and inhomogeneous equation: 
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Here, A are the eigenvalues of the Eq. {TJ, = Wp[2(l — 
j)]^/'' is the dc bias dependent frequency of oscillations 
on the bottom of potential well, [/„((/?). The crossover 
temperature T* is determined by the condition that there 
is the eigenvalue A = 0, 0|. In a global coupling 
case (Eq. Q) the crossover temperature is obtained as 
a solution of the particular transcendent equation: 



-r'-{l-2/N)ujl] = 



(8) 



where ^* = '^^^^'^ ■ Thus, one can see that the crossover 
temperature T* is strongly suppressed for short arrays 
{N ~ 1) for both cases, namely, "inductive" {ojr ^ ujq, 
and a ~ 1) or "capacitive" {ajR ^ loq, and a >> 1) 
types of an external impedance. However, T* recovers to 
the value T* — hujQ/{2TrkB) for long arrays {N » 1). 
Typical dependencies of T*{N) are shown in Fig. 2. 

Now we turn to the MQT regime, where the extremum 
point of the action Seffifn} is the " tau-dependent" in- 
stanton (bounce) solution. At zero temperature and 
in the presence of an external impedance Z, a spatial- 
temporal instanton solution satisfies the equation: 



1 •• / ^ 

—(Pn[r)+a 



dXJ 

dTiGo{T-Ti)iprn~-r^ = 0, (9) 

dip 



where Gq{t) = exp{—u!ji\T\). The solution of Eq.© 
has a following form: a large bounce solution localized on 
a particular junction Z, fiir) = /(r), and a small spatial- 
temporal tail solution distributed over a whole array (see 
schematic in Fig. 1, solid line). The Fourier transform 
of the instanton tail is obtained as 



</5n(w) = - 



(10) 



FIG. 2: Typical dependencies of the crossover temperature 
T*{N) on a number of junctions A'^. Both cases of inductive 
(upper curve) and capacitative (lower curve) impedance for 
particular sets of parameters are shown. 



where go(^) and /(w) arc the Fourier-transform of G'o(t) 
and /(r), respectively. The bounce solution /(r) local- 
ized on the junction I is determined self-consistently from 
the equation: 



'■^p Jo 
Gi{r) 



dnGiiT - ri)/(Ti) -(l-j) + ll^O, 



duj 



(^2 + u;g)e^'^^ 



(11) 



In the absence of a global coupling the instanton solu- 
tion is strongly localized on a particular junction (see 
schematic in Fig. 1, dashed line), i.e Q 



t . f r U 3v/2(l-j) ^ 
cosh [ujot/2) 



(12) 



Substituting (|10|l in the expression ^ for the effective 
action Sfjfand using a perturbative approach (similarly 
toRefs. [ll[l3), i.e. /(r) ~ /o(r), we obtain the MQT 
escape rate (in physical units) as 



where 
X = ■ 

where 



Foexp 



72Ej 



2^/\l-jf/\l + x) 



(13) 



dx 



x'*(x2 + l)gQ(x) sinh '^{nx) 
v^ + l + ^{N-l)x^goix) 



go{x) = l/[x^ -I- (ujn/ujof 



(14) 



(15) 



Here, the parameter Fq is just proportional to N. A 
parameter x having a positive value, characterizes a sup- 
pression of the MQT due to the presence of the charge 
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interaction between Josephson junctions of the array and 
an external shunting impedance. For short Josephson 
junction array {N ~ 1), such a MQT suppression can 
be rather large for moderate values of a. However, as 
we turn to large Josephson junction arrays {N ^ 1) a 
standard MQT behaviour is recovered. Quantitatively 
an enhancement of MQT depends strongly on parame- 
ters a and w^,. The expression ()14(l can be simplified 
in two limits: cjj^ 3> o^o (" capacitative impedance") and 
LUf(, ^ ujQ ("inductive impedance") as 



LUR < UJQ. 



(16) 



Typical dependencies of the MQT escape rate F mqt 
on the dc bias current / for various values of N are pre- 
sented in Fig. 3. One can see a giant enhancement of the 
MQT escape rate as we turn from short to long Joseph- 
son junction arrays. This enhancement results from a de- 
crease of the slope of the bias current dependence escape 
rate. Comparing our theoretical predictions with the ex- 
perimental curves published in Ref. (see Fig. 5 in 
Ref. 0) we find a good agreement for both the crossover 
temperature T* and the dependence of Tmqt {I) for the 
inductive type of a shunting impedance. Therefore, in or- 
der to choose between two models, i.e. a nearest-neighbor 
intrinsic charge interaction or external global charge cou- 
pling, one needs additional independent measurements of 
Debye screening length or to tune the MQT by vari- 
ation of Z. 

In conclusion we have shown that the dissipative (de- 
coherence) effects can be strongly suppressed in long 
{N >> 1) ) Josephson junction series arrays with a global 
charge interaction. The both dissipation and global 
charge interaction can be introduced through an exter- 
nal shunting impedance. This effect manifests itself as 
a giant enhancement of the MQT escape rate for the 
switching from the superconductive state to a resistive 
one (see Fig. 3). A giant MQT enhancement is explained 
through an excitation of spatial-temporal charge instan- 
ton distributed over a whole array. 
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FIG. 3: The dependence of the MQT escape rate Vmqt on 
the dc bias current / for various values of A'^ = 1,50,100. 
Both cases of capacitative (a) and inductive (b) impedance 
for particular sets of parameters are shown. 
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